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Abstract. We consider the centrifugal acceleration of 
charged test particles by rotating magnetospheres which 
are widely believed to be responsible for the relativistic jet 
phenomenon in active galactic nuclei ( AGN) . Based on an 
analysis of forces the equation for the radial accelerated 
motion is derived and an analytical solution presented 
under the assumption of an idealized spherical magneto- 
sphere. We show that the rotational energy gain of charged 
particles moving outwards along rotating magnetic field 
lines is limited in general by (i) inverse-Compton losses 
in the radiation field of the disk in which the magneto- 
sphere is anchored and (ii) the breakdown of the bead- 
on-the-wire approximation which occurs in the vicinity of 
the light cylinder. The corresponding maximum Lorentz 
factor for electrons is of the order of a few hundred for 
the sub-Eddington conditions regarded to be typical for 
BL Lacs. In AGN with enhanced accretion rate the accel- 
eration mechanism seems to be almost inefficient due to 
increasing inverse-Compton losses. 

Key words: acceleration of particles - radiation mecha- 
nisms: nonthermal - gamma rays: theory - galaxies: active 



1. Introduction 

The origin of the nonthermal, highly variable emission 
in active galactic nuclei (AGN) has been widely dis- 
cussed. Several acceleration mechanisms have been pro- 
posed which may explain the observed high energy emis- 
sion extending up to TcV energies at least in three blazars 
(Mkn 421, Mkn 501, lES 2344-^514: e.g. Catanese 1999). 
Fermi-type particle acceleration mechanisms in relativistic 
jets seem to be very effective but require a seed popula- 
tion of electrons with Lorentz factors of at least 100. Up 
to now, it remains a problem to be solved, how this pre- 
acceleration is achieved (e.g. Kirk, Melrose, Priest 1994). 

Since the pioneering work of Gold in the late 1960s 
(Gold 1968, 1969), centrifugal driven outflow of matter 
has often been discussed in the context of pulsar emis- 
sion theory (for recent contribution see e.g. Machabeli 
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& Rogava 1994; Chedia et al. 1996; Gangadhara 1996; 
Contopoulos et al. 1999). In the case of accreting black 
hole systems (e.g. AGN) Blandford and Payne (1982) first 
pointed out that centrifugal driven outflows (jets) from ac- 
cretion disks are possible, if the poloidal field direction is 
inclined at an angle less than 60° to the radial direction. 
In such models, a rotating magnetosphere could emerge 
from an accretion disk or the rotating black hole itself 
(Blandford & Znajek 1977) initiating a plasma outflow 
with initially spherical shape until the flow is coUimated 
on a scale of less than a few hundred Schwarzschild radii 
(e.g. Camenzind 1995, 1996; Fendt 1997). For a rapidly 
rotating black hole system, the critical angle mentioned 
above could be as large as 90° (Cao 1997). 

In magnetohydrodynamical scenarios for the origin of 
relativistic jets, centrifugal acceleration is rather limited, 
leading to maximum bulk Lorentz factors of the order of 10 
(Camenzind 1989). Nevertheless, it seems quite interest- 
ing whether supra-thermal test particles (e.g. from mag- 
netic flares on the accretion disk) could be accelerated to 
even higher energies by such rotating magnetospheres. Re- 
cently, Gangadhara & Lesch (1997) proposed a model for 
spinning active galactic nuclei in which charged test parti- 
cles are accelerated to very high energies by the centrifugal 
force while moving along rotating magnetic field lines. Ac- 
cording to their calculations, the nonthermal X-ray and 7- 
ray emission in AGN could arise via the inverse-Compton 
scattering of UV-photons by centrifugal accelerated elec- 
trons. 

In this paper, we reinvestigate the acceleration of 
charged test particles in an idealized two-dimensional 
model where the magnetic field rotates rigidly with a 
fraction of the rotational velocity of the black hole (cf. 
Fendt 1997). Centrifugal acceleration occurs as a conse- 
quence of the bead-on-the-wire motion. A charged particle 
gains rotational energy as long as it is directed outwards 
but we show that its energy gain is substantially limited 
not only by inverse-Compton losses but also by the effects 
of the relativistic Coriolis force. 

Based on an analysis of forces, the special relativistic 
equation of motion is derived in Sect. 2. This equation 
is solved in closed form in Sect. 3. Sect. 4 gives an esti- 
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mate for the maximum Lorentz factor attainable in the 
case of AGN. The results are discussed in the context of 
the particle acceleration problem for rotating AGN jets in 
Sect. 5. 

2. Analysis of forces in a rotating reference frame 

Usually the motion of a particle along rotating magnetic 
field lines is treated in the bead-on-the-wire approxima- 
tion where a bead is assumed to follow the rotating field 
line and experiences centrifugal acceleration (or decelera- 
tion) while moving in the outward direction (e.g. Macha- 
beh & Rogava 1994; Chedia et al. 1996; Cao 1997). This 
simple approach yields interesting results, though, as we 
will show further below, such an approximation breaks 
down in the region near the light cylinder. 

Let us consider the forces acting on a particle in a ro- 
tating frame of reference (Gangadhara 1996; Gangadhara 
& Lesch 1997). A particle with rest mass m and charge 
g, which is injected at time to and position tq with initial 
velocity vq parallel to the magnetic field line B^ito) ex- 
periences a centrifugal force in the radial direction given 

by 



' cf 



m7 (Jl X r) X rj , 



(1) 



where 7 is the Lorentz factor of the particle and Q, — ^.e^ 
is the angular velocity of the field. Additionally, there is 
also a relativistic Coriolis force in the nonincrtial frame 
governed by the equation 



(2) 



which acts as a deviation-force in the azimuthal direction. 
In the inertial rest frame the particle sees the field line 
bending off from its initial injection position. Hence, it 
experiences a Lorentz force, which may be written as 



-Fl = 9 (^'rcl X B) , 



(3) 



where Urci is the relative velocity between the particle and 
the magnetic field line and where the convention c = 1 is 
used. Due to the Lorentz force a charged particle tries to 
gyrate around the magnetic field line. Initially, the direc- 
tion of the Lorentz force is perpendicular to the direction 
of the Coriolis force, but as a particle gyrates, it changes 
direction and eventually becomes antiparallel to the Cori- 
olis force. Hence one expects that the bead-on-the-wire 
approximation holds, if the Lorentz force is not balanced 
by the Coriolis force. In this case, the accelerated motion 
of the particle's guiding center due to the centrifugal force 
is given by 



7 



dV 
dt2 



dr d7 
dt dt 



(4) 



where r denotes the radial coordinate and 7 — 
l/Vl — ri^r^ — f2. The bead-on-the-wire motion for the 



guiding center breaks down, if the Coriolis force exceeds 
the Lorentz force, i.e. if the following inequality, given by 
the azimuthal components of the forces, holds: 



d7 I f Bq Wrci „ dr 

— > — 2 7 — 

At r \ mfl dt 



(5) 



3. Analytic solution for the radial acceleration 

The general solution of Eq. (^) can be found using the sim- 
ple argument that the energy E of the particle in the rotat- 
ing reference frame is constant. If Eq denotes the energy of 
the particle in the inertial rest frame, then the energy E in 
the uniformly rotating frame (angular velocity Q) is given 
in the non-rclativistic case by: E = EQ — mVl?r^ (e.g. Lan- 
dau & Lifshitz 1960). The generalisation of this equation 
to the relativistic case is straightforward and leads to the 
transformation 



(6) 



where 7 is the Lorentz factor defined above. 

Assume now, that in the general case a particle is in- 
jected at time t = to and position r = with initial 
velocity v = vq. Then, from Eq. (||) it follows that the 
time-derivative of the radial coordinate r may be written 
as 



Ar{t) 
dt 



V(l-17V2)[l-m(l-172r2 



(7) 



where m = {1 — Vt^rl — vl) / {1 — Vt^rlY . In the case tq = 0, 
this expression reduces to the equation given in Henriksen 
& Rayburn (1971). 

The equation for the radial velocity = dr/dt, Eq. (^, 
can be solved analytically (Machabeli & Rogava 1994), 
yielding 



r{t) = -cn(Ao-l^t), 



(8) 



where cn is the Jacobian elliptic cosine (Abramowitz & 
Stegun 1965, p.569ff), Aq a Legendre elliptic integral of 
the first kind: 



Ao 



00 



d0 



(1 - m sin 61)1/2 



(9) 



and where 0o is defined by ipQ — arccos(r2 rg). By using 
Eq. (H), the time-derivative of r can be expressed as r = 
dn(Ao — fit) sn(Ao — Jit). Note that the Jacobian elliptic 
functions sn and dn are usually defined by the relations 
sn^ -I- cn^ — 1 and m sn^ + dn^ = 1 . 

Using Eq. (^, the Lorentz factor may be written as a 
function of the radial coordinate r: 

7 = ^ . (10) 



(11) 



Thus, in terms of Jacobian elliptic functions one gets 
1 

7 = 



m[sn(Ao - r^t)]^ 
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For the particular conditions where the injection of a 
test particle is described by r{to = 0) = and v{to — 
0) = t'o, the time-dependence of the radial coordinate is 
given by a much simpler expression: In this situation, Aq 
reduces to a complete elliptic integral of the first kind 
and therefore, after a change of the arguments, the time- 
dependence of the radial coordinate becomes 



vq sn(f2 1) 
Sldn(r2t) 



If one considers non-relativistic motions, 
and the special case ro = 0, Eq. (||) 
Abramowitz & Stegun 1965) 



where m c 
reduces to 



r{t) = Vq tanh(f7 t)/n sech(f2 1) = vo sinh(il t)/il . 



(12) 

~ 1, 
(cf. 

(13) 



This expression is known to be the general solution of the 
equation 



(14) 



which describes the motion of a particle due to the cen- 
trifugal force in the non-relativistic limit. In Fig. |l|, we 
compute the time-dependence of the radial coordinate r 
for different initial conditions under the (unphysical) as- 
sumption that the bead-on-the-wire motion continues un- 
til the light cylinder (with radius tl) is reached. Note that 
in the relativistic case all particles would turn back at the 
light cylinder due to the reversal of the centrifugal accel- 
eration (e.g. Machabeh & Rogava 1994). 

Using the definition of the Lorentz factor, the equation 
for the accelerated motion, Eq. (|^), may also be written 
as 



dV 



1 - r2 



(15) 



(cf. Chedia et al. 1996; Kahniashvili et al. 1997). By insert- 
ing the above relations, the solution for the radial acceler- 
ation could be expressed in terms of the Jacobian elliptic 
functions: 



¥ = n ■ cn(Ao - r^t) [1 - 2dn^(Ao -Vti)]. 



(16) 



According to our simple model, one expects that a charged 
test particle gains energy due to rotational motion as 
long as it is directed outwards. Therefore the relativis- 
tic Lorentz factor increases with distance r as the particle 
approaches the light cylinder. This is illustrated in Fig. || 
where we plot the evolution of the relativistic Lorentz fac- 
tor 7 as a function of r for different initial velocities wq 
and fixed rg = tl/IO (using a typical light cylinder radius 
of tl — lO^^cm). Note that 7(r/rL) is not scale-invariant 
with respect to the injection velocity vq (i.e. the injec- 
tion energy). If one identifies Eq. ( p^ ) with the general 
expression for the centrifugal force, which reduces in the 
non-relativistic limit to the well-known classical expres- 
sion, the centrifugal force changes its signs and becomes 




Fig. 1. The time-dependence of the radial coordinate r, plot- 
ted for the initial conditions vo = 0.99 c and ro =0.1 tl (solid 
line), Vq = 0.6c and ro = 0.4r_L (short dashed), vq — 0.4c 
and ro — 0.9 (dotted-short dashed); also indicated is the 
non-relativistic limit: r{t) Q/vo = sinh(f2t) (dotted). 



7„=i qq 



7„=iq. 



Fig. 2. The relativistic Lorentz factor 7 for a particle ap- 
proaching the light cylinder using ro — tl/IO and injection 
Lorentz factors 70 — 10 (solid line), 70 = 50 (dotted) and 
70 = 100 (dashed). 



negative for r^/rf^ > 1 — (2m)~^ (see Fig. ^; cf. also 
Machabeli & Rogava 1994). Hence, if one assumes that 
the bead-on-the-wire approximation holds in the vicinity 
of the light cylinder, the radial velocity becomes zero at 
the light cylinder and changes direction in any case. Ac- 



4 



F.M. Rieger & K. Mannheim: Particle acceleration in AGN 



cordingly, a crossing of the light cylinder within the bead- 
on-the-wire approximation, as mentioned for example in 
Gangadhara & Lesch 1997, is not physical (cf. Fig. |l|). 
The reversal of the direction of centrifugal force accord- 
ing to which the centrifugal force may attract rotating 
matter towards the centre is well-known in strong grav- 
itational fields (for Schwarzschild geometry: Abramow- 
icz 1990; Abramowicz & Prasanna 1990; for Kerr geome- 
try: Iyer & Prasanna 1993; Sonego & Massar 1996). For 
illustration, we compute in Fig. the evolution of the 
effective radial acceleration = d r/dt^ as a function of 
the radial coordinate r for different initial velocities. Obvi- 
ously, there exists a point where the effective acceleration, 
i.e. the centrifugal force, becomes negative. 
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Fig. 3. The radial acceleration Or as a function of r/rL for the 
initial conditions ro = and wo = 0.3c (dashed), t^o = 0.6c 
(dotted). 



4. Estimate of the maximum Lorentz factor 

Consider now the acceleration of electrons via rotating 
magnetospheres in AGN. Imagine an electron which moves 
along a rotating magnetic field line towards the light cylin- 
der. Generally, one expects that there are two processes 
which could limit the energy gain of a particle: 

First, there are inverse- Compton energy losses due to 
interaction with accretion disk photons: low energy ac- 
cretion disk photons are scattered to higher energies by 
the accelerated electrons so that the photons gain energy 
while the electrons lose energy. Near the disk the electrons 
might encounter a very strong disk radiation field, which 
substantially limits the maximum attainable energy (this 
need not be the case if electrons are accelerated far away 
from the disk, e.g. Bednarek, Kirk & Mastichiadis 1996). 
The maximum energy, which an electron is able to reach 



under the influence of inverse-Compton scattering is given 
at the point where the acceleration time scale equals the 
cooling time scale. In the case, where the energy of the 
photon in the electron rest frame is small compared to 
the energy of the electron (Thomson scattering), the cool- 
ing time scale for inverse-Compton losses can be approxi- 
mated by (e.g. Rybicki & Lightman 1979) 



tf^,, = 3 X 10^ 



7 



(72-l)C/.. 



ad 



(17) 



where C/rad = t i^Ais^l "^{^ is the energy density of the 
disk radiation field and t < 1- 

If one uses Eq. (|l^), the acceleration time scale tacc may 
be written as: 



2172r^l-m(l-172 ^2) 



(18) 



By equating this two time scales we obtain an estimate 
for the maximum electron Lorentz factor 7max- 

A second, general constraint, which was neither con- 
sidered by Machabeli & Rogava (1997) nor used in the cal- 
culation by Gangadhara & Lesch (1997), is given by the 
breakdown of the bead-on-the-wire approximation which 
occurs in the vicinity of the light cylinder. Beyond this 
point, where the Coriolis force exceeds the Lorentz force 
[see condition Eq. (H)], the particle leaves the magnetic 
field line so that the rotational energy gain ceases. Hence 
the acceleration mechanism becomes ineffective. In the 
case of AGN, where the magnetic field strength is much 
smaller than in pulsars, this constraint may be quite im- 
portant. 

For illustration, we apply our calculations in the fol- 
lowing to a typical AGN using a central black hole mass 
-^BH — 10^ Mq and a light cylinder radius tl — 
lO^^TOg cm, where Mq denotes the solar mass. The Ed- 
dington luminosity, i.e. the maximum luminosity of a 
source of mass Mbh which is powered by spherical accre- 
tion, is given by iEdd — 10'*^ergs s~^. Typically, we may 
express the disk luminosity as Ldisk — le x -^Edd, with 
lO^'' < < 1- Thus, the equipartition magnetic field 
strength at the radius r is given by B{r)^ — 2Ldisk/r'^- 
Electrons are assumed to be injected at an initial posi- 
tion ro ~ 0.4 tl with a characteristic escape velocity from 
the last marginally stable orbit around a black hole of 
vq ~ 0.6 c. By applying the two constraints above, we get 
three generic regimes for the acceleration of electrons by 
rotating magnetospheres: 

1. the region, in which inverse-Compton losses dominate 
entirely over the energy gains, leading to an inefficient 
acceleration (generally in the case of Eddington accre- 
tion, i.e. 4^1). 

2. the region, in which inverse-Compton losses are impor- 
tant but not dominant (generally the sub-Eddington 
range: Ic < 2 x 10"^). In this case the acceleration 
mechanism works, but there exists a maximum Lorentz 



F.M. Rieger & K. Mannheim: Particle acceleration in AGN 



5 



factor given at the position where the energy gain is ex- 
actly balanced by losses. This is illustrated in Fig. ^, 
where we calculate the cooling and the acceleration 
time scale as a function of the Lorentz factor 7 for 
Ze = 5 X 10"'^. For this value, the maximum Lorentz 
factor is roughly 7 ~ 150. Typically, the maximum 
Lorentz factors in this range are of the order of 100 to 
1000 (see Fig. |). 



Lj.^^ = 5x lO''^ ergs 




Fig. 5. Maximum electron Lorentz factor 7max attainable un- 
der the influence of inverse-Compton losses as a function of the 
disk luminosity Ldisk for r = 0.5 (dotted) and r = 1 (solid), 
where r = 4 7rr£l7rad/idisk and Ura.d being the energy density 
of the disk radiation fleld. 



Fig. 4. Cooling times scale t^^^i for inverse-Compton scatter- 



ing, Eq. (|17|), and acceleration time scale tacc, Eq. (^^, as 
a function of the Lorentz factor 7 using Ic = 5 x 10 and 
T — 1. The maximum electron Lorentz factor, given at the 
position where the cooling time scale equals the accelerations 
time scale, is approximately 150. 

3. the region, in which the inverse-Compton losses are 
rather unimportant (generally Ic < 10~^). In this case, 
the maximum Lorentz factor is determined by the 
breakdown of the bead-on-the-wire approximation [see 
Eq. (H)], which yields a general upper limit for the 
Lorentz factor of the order of 1000. This limit is found 
if one approximates Urei by the light velocity which 
amounts to the highest value for the Lorentz forces. 
The results are shown in Fig. ^, where we also allow 
the injection position to vary. We wish to note, that 
the results, presented in Fig. ^, depend essentially on 
the assumed intrinsic magnetic field strength and the 
size of the light cylinder radius (i.e. the angular ve- 
locity). Generally, for a sufficient approximation, the 
maximum Lorentz factor is given by: 

2/3 



1 



7i 



BirL)q 



(19) 



Thus, even if one uses a magnetic field strength of 
-B(rL) = 100 G, which is roughly three times the cor- 
responding equipartition field, the maximum Lorentz 
factor does not exceed 2.5 x 10"^. 



B„-30G 



Fig. 6. Maximum electron Lorentz factor 7max as a function 
of the initial injection position ro for vo = 0.6 c and -B(r'L) ~ 
30 G (i.e. a disk luminosity Ldisk 1.35 x 10*^ ergs s"^). The 
dotted line shows the decrease in efficiency of energy gain e/ = 
7max/7o, while the dashed line indicate the relativistic limit for 
injection given by the condition 1 



rg > 0. 



5. Discussion 

We have considered the acceleration of charged test parti- 
cles via rotating magnetospheres based on a model topol- 
ogy which is motivated by the standard model for AGN 
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(cf. Bcgclman 1994; Camcnzind 1995; Fcndt 1997). Ac- 
cordingly, the jet magnetosphere originates very close to 
the central black hole from an accretion disk, with initially 
spherical profile until the rclativistic jet is collimated to a 
cylindrical shape outside the light cylinder. 

The centrifugal particle acceleration model described 
in this paper extends the calculations by Machabcli & Ro- 
gava (1994) and Gangadhara & Lesch (1997). We find that 
the maximum Lorentz factor attainable for an electron 
moving along a rotating magnetic field line is substantially 
limited not only by radiation losses (e.g. inverse- Compton) 
but also by the breakdown of the bead-on-the-wire approx- 
imation which occurs in the vicinity of the light cylinder. 
Due to these limiting effects, the acceleration of particles 
by rotating magnetospheres seems to be rather less im- 
portant in the case of AGN. Our current calculations in- 
dicate, that for sub-Eddington accreting black holes, such 
as black holes with advection-dominated accretion flows 
(e.g. Narayan & Yi 1994; Narayan 1997), efficient prc- 
acceleration of electrons to Lorentz factors of the order of 
a few hundred might be possible, at least under the highly 
idealized conditions of our analytical toy model. It seems 
interesting that the highest energy gamma rays have been 
discovered from AGN of the BL Lac type which very likely 
accrete in a sub-Eddington mode (e.g. Celotti, Fabian 
& Rees 1998). Under such conditions, inverse-Compton 
scattering of accretions disk photons with energy 0.1 
keV produces gamma rays with a maximum energy of 
~ 100 (7max/10"^)^ MeV, which is in general too low to 
explain the observed high-energy gamma rays in blazars 
(e.g. Kanbach 1996; Catancsc 1999). 

We wish to mention that the results in this paper es- 
sentially depend on the assumed intrinsic magnetic field 
and the angular frequency Q, i.e. on the size of the light 
cylinder radius (7max oc (B/fl)^^^). Therefore, one could 
find a way out of the problem above, for example, by 
assuming a light cylinder radius in BL Lac type objects 
which is much greater than 10^^ cm for a black hole mass 
of Mbh = 10^ Mq. However, in view of magnetohydro- 
dynamic models already existing, this seems to be rather 
improbable. In any case, the acceleration of supra-thermal 
test particles by rotating magnetospheres might possibly 
provide an interesting explanation for the prc-acceleration 
which is required for efficient Fermi-type particle acceler- 
ation at larger scales in radio jets. 

There arc several restrictions on our approach, e.g. we 
have assumed a projected, two-dimensional geometry and 
rigid rotation of magnetic field lines almost up to the light 
cylinder, hence, concerning the last point, neglected a kind 
of toroidal twist (Begelman 1994), when the inertial forces 
overcome the tension in the field line so that the field line is 
swept back opposite to the sense of rotation. However, one 
would not expect that these restrictions alter our conclu- 
sions essentially since they should lower the upper limit for 
the maximum Lorentz factor by making the acceleration 
mechanism ineffective somewhat earlier. Another restric- 



tion is the use of special relativity in our analysis, which 
is only justified far away from the black hole. A detailed 
general relativistic model is needed to assess whether this 
might affect the results very strongly or not. 
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